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■ Using the QCD factorization approach, we reexamine the two-body hadronic charmless 

\o . 

I S-meson decays to final states involving a pseudoscalar (P) and a vector (U) meson, 

^ • with inclusion of the penguin contractions of spectator-scattering amplitudes induced 

Qh! by the b Dg*g* (where D = d or s, and g* denotes an off-shell gluon) transitions, 

• which are of order a^. Their impacts on the CP- aver aged branching ratios and CP- 

- ^ , violating asymmetries are examined. We find that these higher order penguin contraction 

■ contributions have significant impacts on some specific decay modes. Since B irK* , Kp 
decays involve the same electro- weak physics as P — > ttK puzzles, we present a detailed 
analysis of these decays and find that the five R-ratios for B — > vrET*, Kp system are in 
agreement with experimental data except for R{ttK*). Generally, these new contributions 
are found to be important for penguin-dominated B PV decays. 
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1 Introduction 



The study of hadronic charmless 5-meson decays can provide not only an interesting avenue to 
understand the CP violation and the flavor mixing of quark sector in the Standard Model (SM), 
but also powerful means to probe different new physics scenarios beyond the SM. With the 
operation of i?-factory experiments, huge amount of experimental data on hadronic 5-meson 
decays has been analyzed with appreciative precision. To account for the experimental data, 
theorists are urged to gain deep insight into the mechanism of rare hadronic 5-meson decays, 
and to reduce theoretical uncertainties in determining the flavor parameters of the SM from 
experimental measurements. 

In the past years, much progress has been made in understanding the hadronic charmless 
S-meson decays: several novel methods, such as the "naive" factorization (NF) jj, the per- 
turbative QCD method (PQCD) j2j, the QCD factorization (QCDF) jHj, and the soft collinear 
effective theory (SCET) j3], have been proposed; in addition, some model-independent meth- 
ods based on (approximate) flavor symmetries have also been used to analyze the rare hadronic 
5-meson decays |5J E] . These methods usually have quite different understandings of the rare 
hadronic 5- meson decays, and hence the corresponding predictions are also quite different. 
General comparison between these various methods can be found, for example, in Ref. [Zj. 
Since we shall adopt the QCDF approach in this paper, we would only focus on this approach 
below. 

The QCDF approach, put forward by Beneke et al. a few years ago, has been used widely 
to analyze the two-body hadronic i?-meson decays j2l El El UHl HH 1121 HBj • The essence of the 
approach can be summarized as follows: since the h quark mass is much larger than the char- 
acteristic scale of hadronic interaction, Aqcd, to leading power in the heavy quark expansion, 
the hadronic matrix elements relevant to two-body hadronic 5-meson decays can be factorized 
into perturbatively calculable hard scattering kernels and universal non-perturbative parts pa- 
rameterized by the form factors and the meson light cone distribution amplitudes (LCDAs). 
This scheme has incorporated elements of the NF approach (as the leading contribution) and 
the hard-scattering approach (as the sub- leading corrections), and provides a powerful and 
systematical means to compute the radiative (sub-leading nonfactorizable) corrections to the 
NF approximation for the hadronic matrix elements. In particular, the strong phases, which 
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are very important for studying the CP violation in 5-meson decays, are calculable from the 
first principle. Detailed proofs and arguments can be found in Ref. 0, and current status and 
recent developments of this approach have also been reviewed recently in ^3] . 

In a recent work jT3], we have studied the higher order penguin contractions of spectator- 
scattering amplitudes induced by the h Dg*g* transitions (where D = d ot s, depends on the 
specific decay mode, and the off-shell gluons g* are either emitted from the internal quark loops, 
external quark lines, or splitted off the virtual gluon of the penguin diagrams), and investigated 
their impacts on the CP- averaged branching ratios and CP- violating asymmetries oi B ^ 
TTvr, ttK decays. It has been found that these higher order penguin contraction contributions 
are not negligible in two-body hadronic P-meson decays, particularly in the penguin-dominated 
P ttK decays. Thus, combining the findings in the literature fSl QEl ^1 it would 
be worthy to take into account these higher order penguin contraction contributions to the 
exclusive hadronic P-meson decays. This encourages us to further investigate their impacts 
on the hadronic charmless P PV (where P and V denote pseudoscalar and vector mesons, 
respectively) decays. 

P — s> PV decays are closely related to their PP counterparts because of their similar flavor 
structures, however, these modes have apparent advantages in some cases. For example, due 
to the less penguin pollution, P ^ vrp decay modes are more suitable than P ^ tttt ones for 
extracting the weak angle a of the unitarity triangle of the Cabibbo-Kobayshi-Maskawa (CKM) 
matrix [20] • Studies on two-body hadronic P PV decays are therefore very helpful to deepen 
our understandings of the rare hadronic P-meson decays. Earlier theoretical studies on B^PV 
decays based on various approaches can be found, for example, in Refs. 1^ . With the 
accumulation of new experimental data and the theoretical improvements, these B^PV decay 
modes have also been reanalyzed recently [OlEIl^BEl- In this paper, we will reexamine these 
hadronic charmless P PV decays within the framework of QCDF approach and take into 
account the higher order penguin contractions of spectator-scattering amplitudes as mentioned 
above. Here we do not consider the decay modes with an rj or rj' meson in the final states, since 
in this case there are many additional unknown parameters pertaining to these two particles, 
such as their contents, mixing angles, and the anomaly g — g~r]^'^ coupling, which would hinder 
us from getting reliable theoretical predictions. 
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This paper is organized as follows. Sec. |2 is devoted to the theoretical framework. In 
this section, we first give the relevant formulas describing the decay amplitudes of hadronic 
B PV decays at next-to- leading order in as, and then take into account contributions 
of the higher order penguin contractions of spectator-scattering amplitudes induced by the 
b — Dg*g* transitions. In Sec. El we give our numerical results for CP-averaged branching 
ratios and CP- violating asymmetries, and discuss the impacts of the higher order corrections 
on these quantities. Detailed analysis of the interesting decays B ttK* and B Kp, are 
also presented in this section. Finally, we conclude with a summary in Sec. |3] Some useful 
functions and the input parameters used in this paper are collected in Appendix A and B, 
respectively. 



2 Theoretical framework for B PV decays 

2.1 The effective Hamiltonian for hadronic 5-meson decays 

Using the operator product expansion (OPE) and the renormalization group equation (RGE), 
the low energy effective Hamiltonian for hadronic charmless P-meson decays in the SM can be 
written as [231 



T^es — —7= X! {CiQ\ + C2Q2+ X! Gi Qi + Q-j^ -I- C^g Qs 

V p=u,c \ i=3,...,10 

where \p = Vpi, V*^ (for b ^ s transition) and = Vpb V*^ (for b d transition) are products 
of the CKM matrix elements, and the unitarity relation — Af^ = A^'-* -|- A^'^ has been used. The 
effective operators, Qi, governing a given decay process, can be expressed explicitly as follows. 

(i) Current-current operators: 

Qi = ipb)v-A{Dp)v-A , Q2 = iPibj)v-A{DjPi)v-A , (2) 

(ii) QCD-penguin operators: 

Q3 = {Db)v-A {qq)v-A , Qa = {Dibj)v-A {qjqi)v-A , 

Qb = iDb)v^A Y iqq)v+A , Q& = {Dibj)v-A Y iqjqi)v+A , (3) 
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(iii) Electroweak penguin operators: 

Qj = {Db)v~A |eg(gg)y+A , Qs = {Dibj)v-A hqiQjQi)v+A , 

<? <? 

Qg = {Db)v-AY1 Qio = {Dibj)v-Aj2 ^'^qi^jlih-A , (4) 

Q Q 

(iv) Electro- and chromo-magnetic dipole operators: 

Q^i = ^2'^bDa^,{l + 75) F^^'^h , Qsg = ^m,D a,, (1 + 75) G^'^h , (5) 

where {qiq2)v±A = ^i7m(1 ^15)Q2, hj are colour indices, is the quark electric charge in units 
of |e|, and a summation over q = u,d,s,c,b is implied. For b ^ d transition induced decay 
modes, D = d, while for b ^ s transition induced ones, D = s. 

The Wilson coefficients Ci{fi) in Eq. represent all the physics contributions higher than 
the scale /i ~ 0{mb). Numerical results for these coefficients evaluated at different scales can 
be found in Ref. j2Sl- 

2.2 Decay amplitudes at next-to-leading order in ag 

With the low energy effective Hamiltonian given by Eq. (^, the decay amplitude for a general 
hadronic charmless B PV decay can be written as 

{pv\n^^\B) = % E {PvmB) . (6) 

V p=u,c 

Then, the most essential theoretical problem in the calculation of the decay amplitude resides in 
the evaluation of the hadronic matrix elements of the local operators contained in the effective 
Hamiltonian, {PV\Q^\B) . With the QCDF approach, they could be simplified to a large extent. 
To leading power in Aqcd/''^;), but to all orders in perturbation theory, these hadronic matrix 
elements obey the following factorization formula [S] 

+ Tf" *fB^B* fp $p * /y $y , (7) 

where $jvf are the LCDAs of the meson M, the * products indicate convolutions of the LCDAs 
and the hard-scattering kernels T/'^^ F^'~*^ and denote the heavy-to-light B —> P and 



Figure 1: Graphical representation of the factorization formula. Only one of the two 
form-factor terms in Eq. (|7j) is shown for simplicity. 

iiij Ml M2 M2 




(<:) (t) (g) (li) 

Figure 2: Order corrections to the hard-scattering kernels Tl^^ (coming from the 
diagrams (a)-(f)) and T/^ (coming from the last two diagrams). 

B ^ V transition form factors, respectively. A graphical representation of this formula is 
shown in Fig. d 

When the power suppressed C(Aqcd/"^;)) terms are neglected, T/'^^ are dominated by hard 
gluon exchanges, and hence calculable order by order in perturbative QCD. The relevant Feyn- 
man diagrams contributing to these hard-scattering kernels at next-to-leading order in as are 
shown in Fig. |21 The kernel ^ starts at tree level and, at next-to-leading order in a^, con- 
tains the sub-leading "nonfactorizable" corrections coming from the vertex-correction diagrams 
Figs. Efa-d) and the penguin diagrams Figs. Efe-f). The kernel T/^ contains the hard "nonfac- 
torizable" interactions between the spectator quark and the emitted meson M2. Its lowest order 
contributions are of order and can be depicted by the hard spectator-scattering diagrams 
Figs. 121^ g-h). At leading order, T^^j = 1, T/^ = 0, the QCDF formula reproduce the NF results. 

As stressed in Ref. jH], it should be borne in mind that the factorization formula, Eq. (jZj), 
does not imply that the hadronic 5-meson decays are perturbative in nature. Dominant soft 
contributions to the decay amplitude do exist. However, all these nonperturbative effects either 
are power-suppressed by Aqcd / f^b or can be factorized into the transition form factors and the 
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meson LCDAs. 

With the above discussions about the effective Hamiltonian for hadronic 5-meson decays 
and the QCDF formula for the hadronic matrix element, the decay amplitude for a general 
hadronic charmless B PV decay, in the heavy quark limit, can then be rewritten as 



G 



10 



F ' 



(8) 



p=u,c i=l 

where {PV\Qi\B) p is the factorized hadronic matrix element, which has the same definition as 
that in the NF approach. All the "nonfactorizable" effects are encoded in the coefficients af , 
which are process dependent and can be calculated perturbatively. Following Beneke et al. |12j . 
the general form of the coefficients (z = 1, 10) at next-to- leading order in a^, with Mi 
being the meson picking up the spectator quark and M2 the emitted meson, can be written as 



Cj±lN,r/ji^\ , Ci±lCFOis 



-)Ni{M2) + 



+ PHM2) , (9) 



where the upper (lower) signs apply when i is odd (even). The quantities Vi{M2) account for 
one-loop vertex corrections, Hi{MiM2) for hard-spectator interactions, and Pf (M2) for penguin 
contributions. Explicit expressions for these quantities can be found in Ref. fT^ . 

It is noted that, in calculations of the decay amplitudes for hadronic charmless i?-meson 
decays, the coefficients {i = 3, 10) always appear in pairs. So, for the two-body hadronic 
charmless B — > PV decays, one can define the following quantities af in terms of the coefficients 
af defined in Eq. © [121 

ai(MiM2), 

a2(MiM2), 

(M1M2) - (M1M2) ; 



ai(MiM2) 
a2{MiM2) 

«I(MiM2) 
af(MiM2) 
al,jMiM2) 
al,jMiM2) 



af(MiM2) + af(MiM2); 

<(MiM2) + rf 2 ag(MiM2) 
al{MiM2) - rf 2 a^(MiM2) 

(M1M2) - a?(MiM2) ; 



if M1M2 = VP, 
if M1M2 = PV, 

i iiMiM2 = PV, 
; ii MiM2 = VP, 

if Ml Ms = VP, 



(10) 



a^(MiM2) + a?(MiM2) ; if M1M2 = PV . 



a?o(MiM2)+rf2ag(MiM2); 
a?o(MiM2)-rf ag(MiM2); 



if M1M2 = PV, 
if M1M2 = VP , 



Figure 3: The weak annihilation diagrams of order as- 
with the scale-dependent ratio r^"^ defined as 

where all quark masses are running current masses defined in the MS scheme, and fy{fi) 
is the scale-dependent transverse decay constant of vector meson. Although all these terms 
proportional to r^^ are formally power suppressed by AqcD/'^b in the heavy-quark limit, they 
are not small numerically. In particular, the factor r^{n) is chirally enhanced and important 
for charmless B decays jHl • 

According to the arguments in [Hj, the weak annihilation contributions to the decay ampli- 
tudes are power suppressed, and hence do not appear in the QCDF formula, Eq. (jZj). Neverthe- 
less, these contributions may be numerically important for realistic i?-meson decays. At order 
0{as), the annihilation kernels arise from the four diagrams shown in Fig. 121 They result in a 
further contribution to the hard scattering kernel T/^ in the QCDF formula, Eq.Q. However, 
within the QCDF formalism, these annihilation topologies violate factorization because of the 
end-point divergence. In this work, following the treatment of Refs. jHl 123 , we will introduce a 
cutoff to parameterize these contributions and express the weak annihilation decay amplitudes 
as 

^-n,(5 ^PV)^^J2T. >^^fB /m. /m. h{M,M^) , (12) 

V ^ p=u,c i 

where and /m are the decay constants of the initial B and the final-state mesons, respectively. 
The parameters hi{MiM2) describe the annihilation contributions and their explicit expressions 
can be found in Refs. jHllEl- 

The explicit expressions of the decay amplitudes for hadronic charmless B PV decays, 
including the weak annihilation contributions, can be found, for example, in Refs. ^3 E] • It 
should be noted that, within the QCDF framework, all the "nonfactorizable" power suppressed 
contributions except for the hard spectator interactions and weak annihilation contributions are 
neglected. In addition, in the evaluation of the hard spectator and weak annihilation terms, the 
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running coupling constant and the Wilson coefficients should be evaluated at an intermediate 
scale Hh ~ (AqcD^^b)^^^ rather than the scale fx ~ m^. However, the evolution of Cj(/i) down to 
Hh is highly nontrivial, since the RGE will change below the scale mi). To deal with this problem, 
one may have to turn to SCET which is the appropriate effective theory for QCD below the 
scale. However, in this paper, we restrict ourself to QCDF and adopt the treatments of 
evolution of the Cj(/i) as done in [H], i.e., we do not take into account the charm and bottom 
threshold and evolve the Wilson coefficients in a 5-flavored theory. With this approximation, 
in the evolution of the Wilson coefficients, all logs of the form log fi/M^ have been summed, 
while logs of the form log fi/rrif, and log fi/rric are not. Since the latter two terms are never large 
with /i > mb/2, the approximation would work to the precision in this paper. ^ Specifically, we 
shall use fit = v^A/TyU with Ah = 0.5 GeV in our numerical calculations. 

2.3 Penguin contractions of spectator-scattering amplitudes and their 
contributions to B ^ PV decays 

At the quark level, the h Dg*g* transitions can occur in many different manners as depicted 
by Figs. EHEl For example, one of the two off-shell gluons can radiate from the external quark 
line, while the other one comes from the chromo-magnetic dipole operator Qsg as Figs. El^b) 
and El^c) or from the internal quark loops of the penguin diagrams as Figs. Efb) andlHfc). 
On the other hand, the two off-shell gluons can also radiate from the internal quark loops as 
Figs, ini^d) andinfe) or split off the virtual gluon of the penguin diagrams as Figs. Efa) andlBfa). 
Here we do not consider the Feynman diagrams of the category shown in Fig. |3J since their 
contributions can be absorbed into the definitions of heavy-to-light transition form factors as 
Figs, m^a) andlH^b) and the meson LCDAs as Figs. IHe), or are further suppressed by as 
Figs, life) andm^d). It is easy to clarify this point by comparing the strength of Fig. IH^c) to 
that of Fig.Ei;a). 

As shown by Figs. El and IHl these Feynman diagrams should be the dominant sources con- 
tributing to the penguin contractions of spectator-scattering amplitudes of order a^, since they 
are not two-loop QCD diagrams, and hence there are no additional suppression factor in 
their contributions compared to the genuine two-loop ones of order a"^. Studying these contribu- 
^We thank M. Bcneke for pointing out this point to us. 
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(d) (e) 

Figure 4: Representative Feynman diagrams induced by the h — > Dg*g* transitions which are 
not needed to evaluate. Only the chromo-magnetic dipole operator Qsg contributions are shown. 
With the operator Qsg replaced by the other operators, the corresponding Feynman diagrams 
can also be obtained. 

M2 J1/2 M2 






(a) (b) (c) 

Figure 5: Chromo-magnetic dipole operator Qsg contributions induced by the b 
transitions. 



Dg*g* 



tions could be very helpful for our understandings of the higher order perturbative corrections 
to the rare hadronic 5-meson decays within the QCDF formalism. 

We start with the calculations of the Feynman diagrams in Fig. El In this case, the b quark 
weak decay is induced by the chromo-magnetic dipole operator Qsg, and the calculation is 
straightforward with the result given by 



A, 



QSg 



ATS 
c 



Jo 



X J^dudv ■!^^M2{u)^Mi{v) 



1 



2{l-u){l-v)v 6{l-u){l-v) 
3{3 — u — V + uv) 2 — u 



2{l-uy{l-v)v Q{l-u)u{l-v) 
1 3{3-v) 



Q{l-u){l-v 
3{3 — u — V — uv) 



+ 



2{l-u){l-v)v 
1 + u 



(13) 



2{l-u){l-v)v 6(l-n)(l-'u) 
when Ml is a pseudoscalar and M2 a vector meson. For the opposite case of a vector Mi and 
a pseudoscalar M2, one needs only change the signs of the last two terms in the bracket of 
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Figure 6: Penguin operator Qi contributions induced by the b Dg*g* transitions. 

Eq. ((131). Here At = VtbV*^ and = VtbV*^ are products of the CKM matrix elements, $a/ and 
$m denote the leading-twist and twist-3 LCDAs of the meson M in the final states, respectively. 
In our calculation, we use the LCDAs in the asymptotic limit 



$p(x) = = 6x(l - x). 



1, $^ = 3(2x-l^ 



(14) 



and have neglected the tree-particle LCDAs and deviations from the asymptotic limit. 

In calculating the Feynman diagrams in Fig. IHl we adopt the method proposed by Greub 
and Liniger ^7]. We first calculate the fermion loops in these individual Feynman diagrams, 
and then insert these building blocks into the entire Feynman diagrams to obtain the final 
results. In evaluating the internal quark loop diagrams, we shall adopt the naive dimensional 
regularization (NDR) scheme and the modified minimal subtraction (MS) scheme. In addition, 
we shall adopt Feynman gauge for the gluon propagator throughout this paper. The gauge 
invariance will be guaranteed when the full set of Feynman diagrams are summed with the 
external quarks on- mass-shell However, we must be care of the gauge dependence in 

our calculation, since only a subset 0{a1) Feynman diagrams are calculated. After careful 
checking, we find that each Feynman diagram in Fig. 5. and 6 is gauge independent. The detail 
checking can be found in Appendix C. Analogous to the calculation of the penguin diagrams 
in Fig. I21^e), we should also take into account the two distinct penguin contractions of the 
four-quark operators in the weak interaction vertex . 

As shown in Fig. IHl the first three Feynman diagrams have the same building block /^(/c) 
(corresponding to contractions of the operators Qi,3,9) or /^(/c) (corresponding to contractions 
of the operators Q4,6,8,io)- These building blocks can be depicted by Fig. [7|and given by 



in- 



r(-) (2 - e) (4vr/i2)f ^k, ^ - Pj,) (1 - 75) 
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s,d 



Figure 7: Building bloclis I?^{k) (corresponding to contractions of the operators (5i,3,9) aud 



(corresponding to contractions of the operators (54,6,8, lo) for Figs. |ni^a)-inic). 



X 



dx ■ 



X (1 — x) 



irt^ — x(l — x)kP' — i5 



(15) 



^F(^)(4vr/i2)t(A;,^-A:^7.) 

x I dx 
Jo 



75) T'^ 



x(l 



x] 



m'i — x{l — x) k"^ — i6 



(16) 



where k is the momentum of the off-shell gluon, T"" = with A" the Gell-Mann matrices, gs 
is the strong coupling constant, and nig the pole mass of the quark propagating in the fermion 
loops. We have used d = A — e. After performing the subtraction with the MS scheme, we get 



9s 



^ - ^ In — - G{sg,r) 
3 3 nib 

In — - G{sg,r) 
3 nih 



(k, ^ - k'^,) (1 - 75) 
ik,^-e^^)il-^,)T\ 



(17) 
(18) 



with the function G{sq,r) defined by 

G{sq,r) = — 4 / dxx{l — x) \n[sq — x(l — x)r — i6], 

J 



(19) 



where Sg = m^/m^, r = k'^/ml, and the term i6 is the "e-prescription" . The free indices fi and 
a should be contracted with the gluon propagator when inserting these building blocks into the 
entire Feynman diagrams. 

The sum of the fermion loops in the last two diagrams in Fig. IHl are denoted by the build- 
ing block J°,l{k,p) (corresponding to contractions of the operators (51,3,9) or J^l{k,p) (corre- 



sponding to contractions of the operators (54,6,8,10)5 which are depicted by Fig. |H1 Using the 
decomposition advocated by fOlEI, these building blocks can be expressed as 



T;^^{k,p) {T^,T'} + T;^{k,p) 
f;^{k,p) {^^^^} + f;,(A:,p) 



(20) 
(21) 
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Q4.6,S, 





g'{lJ.,a,k) 



Figure 8: Building blocks J'^{k,p) (corresponding to contractions of the operators (51,3,9) and 
Jutik,p) (corresponding to contractions of the operators <54,6,8,io) for Figs.lHtd) andlH^e). 



where the first (second) part is symmetric (antisymmetric) with respect to the color structures of 
the two off-shell gluons. Here k (p), a (b), and fi (z/) are the momentum, color, and polarization 
of the off-shell gluons, respectively. Below we refer to the gluon with indices (z/, b, p) as the one 
connected to the spectator quark. 

In the NDR scheme, after performing the (shifted) loop momentum integration, we can 
represent the quantities T^^(/c,p) and T^j,(/c,p) as p!!)| IT7j 

k,. 



T^^k.p) 



T~u{k,p) 



Air 

~E{fi, k, p) Ai24 - E{i^, k, p) -p^ Az25 - E{u, k, p) 7^ Az26 



E(/i, z/, k) Ais + E(/x, I/, p) Aid - E{^, k, p) Ai23 

k ■ p 



as 
Air 



k ■ p k ■ p k ■ p 

^ g^y Az2 + i g^,u Azs + 7^ A;^ Azg + 7m Pu Azg + 7^ k^, Azn 



L, (22) 



-7^p^ A«i2 + It^— 37 A«i5 + |(;^ I\i ' ^" ' " 



h h 

A« 

k ■ p 



T'{k,p) + 



19 



A; ■ p 



k ■ p 



16 



k ■ p 



Ai 



18 



Air 



k ■ p 



Ai2o +^7^Ai2i +^^^Ai22 



l^^ky 



k ■ p 



k ■ p 

L, (23) 
(24) 



luk^ 



(25) 
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where L = 1 — 75, and the matrix E in Eq. ()22j) is defined by 

E{fi,u,k) = 'j^'jult -j^K +-f^k^ - Itg^^ 



(26) 



with the second fine obtained in a four dimensional context using the Bjorken-Drell conventions. 
The parameter a in Eq. (j^^ denotes the chiral structure of the four-quark operators in the 
weak vertex with a = ± corresponding to {V — A)®{y ^ A), respectively. Explicit expressions 
for the dimensionally regularized Ai functions can be found in Appendix B of Ref . ^H] ■ 

Equipped with these building blocks, we can now evaluate all the Feynman diagrams in 
Fig. El After direct calculations, the final results of these penguin contractions of spectator- 
scattering amplitudes for hadronic charmless B PV decays can be expressed as 



Qs 



An., 



Ac 



2 4 u 

- - - In G{sp,u) 

3 3 mfo 



fiiu,v) 



+ 



- - - In G{sp,uv) 

3 3 nib 



f2{u,v) + fi{u,v,mp) I , 



(27) 



c 

inA_G(0,^)-G(l,^) 
3 3 rrih 

- - - - GiQ.uv) - G{l,uv) 

3 3 

2 ) 



f2{u,v) 



(28) 
(29) 



^ \\ ^ ^ I dudv I [(n, - 2 



/2(u,w, 0) + f2{u,v,mc) 



+ 



In - (n. - 2) G(0, u) - Gis^, u) - G(l, u) 
3 rUb 



^ In - (uf -2) G(0,uv) -G(sc,uv) -G(l,uv) 
3 nib 



f2{u,v) 



fi{u,v) 
(30) 



Aq, = I "'^^^^^^"'^ \V l^'dudvl^linf - 2) Mu, V, 0) + Uu, v, mj 



+f3{u,v,mb)] + 



4n/ , fi 



In — - (nf - 2) G(0, u) - G(s,, u) - G(l, u) 
3 nib 



+ 



-i^ln^ 
3 nib 



{nj -2)G{0,uv) -G{sc,uv) -G{l,uv) 



f2{u,v) 



fiiu,v) 
(31) 



• Oil Ib fhh /m2 Ai) 







+ 



1 



dudv 



f3{u,v,m^) - -fs{u,v,mb) 



1 In A -Gis,,u) +^Gil,u) 
3 TUb 2 
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+ 



"^In— G{sc,uv) + ^G{l,uv) 

6 rrii, 2 



f2{u,v) 



A 



Qio 



Jo 



e Jo 



+ 
+ 



1 In A -G{s,,u) +Ig{1,u) 
6 mi, 2 



In— G{sc,uv) + ^G{l,uv) 

6 mi, 2 



dudv 

fiiu,v) 



f2{u,v,m^) - -f2{u,v,mi,) 



(32) 



(33) 



with the subscript Qi denoting the contraction of Qi operator in the weak vertex, and 



fi{u,v) 



fi{u,v,mg 



f2,{3){u,v,mg) 



2U + V-3 3{2u + v-3) 

\2 "I 



+ rfi$M,(n)$^,(i; 



12{l-u){l-vy A{l-u){l-v)v 
3{v-3) 1 



A{l-u){l-v)v 12{l-u){l-v) 



2u-l ^3{2u + v -2uv -3j 



12{l-u){l-v) A{l-u){l-v)v 
v-3 3{v-3) ' 



12{l-vy A{l-v)v 



(34) 



f2{u,v) = <i>M,{u)<^M,{v] 



— -2 - rf 1 $M.(n) ^mAv) 

6 (1 — V) 



1 



6 n (1 — v) 



+rf rf$.,(n) $.,(.) 



rf $„,(n)<|.M,(^)-— ^ , (35) 
6{1 — V) 

3Ais ^ 7AiQ 



1 



3At2 



{l-u){l-v) 8{l-u)v 2A{l-u)v 



^ 3Ai8 ^ 7Az23 ^ 7(l-n + t;)Az5 



8{l-v)v 2A{l-v)v 2A{l-u){l-v)v 



3 (Aia - Ais + Ain] 

12{l-u){l-v)^ 8(1-m)(1-w) 



^ 7 (Aze + At2e) ^ 3 (A zg + 2 Ati2 + Azai) 



24(1 -u)v 



+ rf rf*^, («)$„, (^;) 



i(l -n)i; 

7 Ais ^ 3(Ai 2 + Aig - Aiu + Aii7) 



7uAi23 



12(l-t;) ' 8?; 
3 {Ai2 + Ais + Ain) 3u{Ai3 + Ai2i) 



12(1 -m) {1~v) 
7 {Aid + Aise) 



(l-u){l-v) 



{1 — u) V 



24(1 - n)z; 
rf <l>^,(w) <l>M,(t;^ 



3(Az2-Az8) 7(Az23 + 2Ai5) 



8{l-v)v 
{3-2u-2v) 



+ 



24:{l-v)v 
3Ai2 



(36) 



3A23 



2{l-u){l-v)v 8(l-M)(l-t;) 8{l-u)v 
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7Ai 



6 



3 Azg 



24 (1 -u)v t 



± 



7Ai 



23 



± 



7(1 -u + t;) A?5 
(1 - i;) ^ 24 (1 - w ^ 24 (1 - u) (1 - t;) w 

3(A23 + 2Azi2 + Ai2i^ 



+ 



2A(l-u)v 



± 



2(1 -m) (1 ' 8(1 

7Ai5 3 (Ai2 - Aig + Aiir) 



12(l-n)(l-i;) 



(l-n)(l-^) 



mi ' 



2n-2t; + 2nw 



7 M Ai23 



2(l-n)(l-tOf 12(l-u)(l-i;) 



7A25 



3m(A23 + A22l) ^ 7(A26 + A^26) 



12(l-t;) 8{l~u)v 24(1-M)t; 
3 (Ai2 + Ais + Ain) 3 (Ai2 + Ais - Ai^ + Ain) 



^ 8(1-m)(1- 
- rf $„,(n)<|.M,(^;) 



3(Az2- Azg) ^ 7(Az23 + 2Az5) 



(,37) 



2{l-v)v 8{l-v)v 24(l-t;)t; 
when Ml is a pseudoscalar and M2 a vector meson. For the opposite case, i.e., Mi is a vector 
and M2 a pseudoscalar meson, one needs only change the signs of the last two terms in the 
functions fi defined above. At this stage, the Aj functions appearing in Eqs. (jHUI) and ()37p are 
the ones that have been performed the Feynman parameter integrals, whose explicit expressions 
can be found in Appendix B of Refs. fSl ^] • For convenience, we also list them in Appendix 
A. 

With the individual operator contribution given above, the total contributions of the pen- 
guin contractions of spectator-scattering amplitudes can be written as 



1 



^ Ci Aq^ + {C3 - - Cg) Aq^ + C4 Aq^ + Cq Aq^ 



jp=u,c 



+Cs Aq^ + Clo ^Qio + ^ 



89 Q8g 



(38) 



where the superscript '/' indicates the one to be distinguished from the next-to- leading order 
results given by Eqs. (jH)) and (jl2|l . The total decay amplitude is then given as 



{PV\ncs\B) = A{B PV) + A^^'^iB PV) + A^{B PV) . 



(39) 



3 Numerical results and discussions 

With the theoretical expressions given above and the input parameters collected in Appendix B, 
we can now evaluate the branching ratios and CP- violating asymmetries for two-body hadronic 
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charmless B PV decays, with P = (vr, K) and V = {p,uj, K*, 0). For each quantity, we first 
give the resuhs at next-to- leading order in a^, and then take into account the higher order 
penguin contractions of spectator-scattering amplitudes induced by the b Dg*g* transitions. 
The combined contributions of these two pieces, denoted by 0{as + a^), are then given in the 
last. For comparison, results based on the NF approximation are also presented. All the exper- 
imental data are taken from the home page of the Heavy Flavor Averaging Group (HFAG) |26j . 

In order to show the renormalization scale dependence of the branching ratios and CP 
asymmetries, we give results of two cases for each decay mode with the first one evaluated at 
the scale fi = rrih, while the second at the scale fi = mf,/2. In addition, our calculations depend 
on many input parameters, which cause quite large theoretical uncertainties. We will consider 
the main theoretical uncertainties arising from the strange-quark mass (with the ratio rriq/ms 
fixed, all chiral enhancement factors depend on this mass), CKM matrix elements, form 
factors, and the first inverse moment of the 5-meson distribution amplitude A^. ^ 

3.1 Numerical analysis of penguin contractions of spectator-scattering 
amplitudes 

Before presenting numerical results for branching ratios and CP asymmetries, we would discuss 
the relative strength of each Feynman diagram shown in Figs. El and IHl For convenience, 
we denote the decay modes with the pseudoscalar meson picking up the spectator quark by 
B PV, while for the vector meson picking up the spectator quark by -B ^ VP. 

Firstly, we study the relative strength of the three Feynman diagrams shown in Fig. El Since 
contributions of these diagrams are all proportional to 



we have factorized 5*1 out off the numerical results shown in Table ^ 

From the numerical results for the dipole operator Qgg contractions given in Table ^ we 
^ Since the theoretical uncertainties coming from the weak annihilation and twist-3 hard-spectator interaction 
contributions (parameterized by the quantities Xh and Xa) are already known to be quite large [T^. we do not 
consider these uncertainties here and simply use the default values given by Xh{A) ~ ^ogms/ Xh as specified in 
Appendix B. Uncertainties coming from the other input parameters are generally small and have been neglected. 




(40) 
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Table 1: Numerical results of each Feynman diagram shown in Fig.^lwith the asymptotic 
forms of the meson LCDAs. Terms involving the twist-three LCD As are given in unit of 
the factor r^' defined by Eq. (fTT|) . The subscripts Mi = P{V) and M2 = V{P) for the 
B PV ( B VP ) rows. The same for mi 2. 







Decay mode 












r.l^a) 


B ^ PV 


-67.50 


-125.76 


-9.64 


-18.94 


Fij 
















B ^VP 


-67.50 


4.82 


34.71 


-3.79 


Fij 


^s. Efb+c) 


B ^ PV 
B ^VP 


-1.50 
-1.50 


-3.54 
-1.61 


-1.07 
1.86 


-0.42 
0.42 



can see that the main contributions come from Fig. Efa), and the other ones play only a minor 
role. It is also noted that these amplitudes do not have strong phases. 

To analyze strong phase sources and the relative strength of the individual Feynman diagram 
shown in Fig. IHl we study the Ql contraction in the weak vertex. The contribution of each 
Feynman diagram is proportional to 

e . «s /b /mi /m2 w/) /"^^f "^f (0 /.IN 
S2=^ ^1 X^^^' (41) 

which has also been factorized out. The numerical results given in Table El are independent of 

From the numerical results given in Table |2l we have the following observations: (i) contri- 
butions of Figs, in^b) and El^c) are generally much smaller than those of the other three ones, 
and the main contributions come from the diagrams Figs. Efd) andEl^e); (ii) although each 
term labelled by the meson LCDAs in each Feynman diagram has a large imaginary part, and 
hence a large strong phase, the total strong phase of each Feynman diagram is small due to 
cancellations among the four terms, (iii) for each term labelled by the same LCDAs, there also 
exist cancellations between the contributions of the diagrams Fig. IHl^a) and Figs. IHl^d-l-e). 

Thus the total strong phase is found to be quite small after summing all the five diagrams 
shown in Fig. IHl Moreover, the cancellation does not depend on the parameters in 82- 
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Table 2: Numerical results of each Feynman diagram shown in Fig. IHl Others are the 
same as Table [H 







modes 














;.ia:a) 


r V 


— i.oy — iz.oo I 


n 1 7 


1/1 1 n 1 


— U.iO -|- iO.OO I 


U.iZ -|- iO.Oi i 


Fk 














VP 


-1.39 - 12.65Z 


-0.02- 


f 1.28z 


-0.12 + ll.llz 


-0.01 - 0.44 z 






PV 


-0.01 - 1.05i 


-0.12 - 


- 1.2H 


-0.62 + 0.8H 


-0.18 - O.lli 




;s. ini^b+c) 


















VP 


-0.01 - 1.05 z 


-0.39- 


- 1.25z 


-0.08 + 0.78 z 


-0.10-0.19z 




IS. Efd+e) 


PV 


-9.03 + 14.94 z 


19.19 4 


- 28.30 z 


4.32 -21.29 z 


10.82 - 15.69 z 


















VP 


-9.03 + 14.94 z 


14.26- 


(-9.04Z 


0.83 - 16.78 z 


-0.39 - 3.46 z 



3.2 Branching ratios of B ^ PV decays 

In the following discussions, we classify the two-body hadronic charmless B PV decays into 
two categories: the strange-conserving {AS = 0) and the strange-changing {AS = 1) processes. 
The higher order penguin contractions of spectator-scattering amplitudes are expected to have 
more significant impacts on the AS = 1 processes than on the AS* = ones, due to the CKM 
factor suppressions in the latter. 

Numerical results of the CP-averaged branching ratios for these decays are collected in 
Tables El m andEl where the theoretical error bars are due to the uncertainties of the input pa- 
rameters CKM elements, quark masses, transition form factors, and as collected in Appendix 
B. Generally, the theoretical uncertainties are quite large, which are larger than 0{a'^) correc- 
tions for tree-dominated decay modes, but comparable for strong penguin-dominated decay 
modes. For most decay modes, the corrections reduce the renormalization scale dependence 
of the theoretical predictions. 

For AS = decays, since the b ^ d penguin amplitudes are suppressed by the CKM factor 
compared to the b s penguin amplitudes, most of them are expected to be dominated 
by the tree amplitudes, however with a few exceptions. From the numerical results given in 
Tables El and lU we have the following general remarks: 

• The decays 5^ — ir^p^ and Tr'^p^, it~ p^, tx~uj. Our results are generally consistent 
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Table 3: CP-averaged branching ratios (in units of 10^^) of tree-dominated B PV 
decays with AS* = 0. B-^ and B^^"" denote the results without and with the annihilation 
contributions, respectively. Results in columns 0{as + «s) are the ones with the higher 
order penguin contraction contributions included. For each decay mode, the first row is 
evaluated at the scale fi = rrih, while the second one at the scale fi = mfe/2. The theoretical 
errors correspond to the uncertainties of the input parameters collected in Appendix B. 
The NF results are also shown for comparison. 



Decay mode NF B^ EXP. 









0{a,) 


0{as + al) 


0{as) 


0{as + al) 




B- - 




8.76t|i 


8.15tlil 


8.O2I2I80 


1 Q+3.53 
0. 10 2, 63 


8-01111 


8.71}:° 






7 C9+3.36 
' — 2.45 


7.45^y? 


7.36l2,67 




7.3611^0 




B- - 




13.91+^i^ 


13.05+^i 


13.31+^-?^ 


1 Q 00+5.94 


13.48+^^^ 


10.8+}-^ 

— 1.0 






13.08^^^^ 


i2.82i^:i 


13.011^:81 


13.0011^^ 


13.20l|^^ 




b"- 




19.78t?:i 


19.37t?;i 


19.73l^°2f 


20.341^092° 


20.721?:!^ 


13.911? 








20.22^^^11° 


20.481^^65^ 


21.251^^2^3 


21.52l^°8f 










10.51111 


10.471^:^° 


ii.iO_3 g2 


ll-lllt?' 


10.11?:^ 










10.86lii 


11.5711^^ 


1 1 c;9+4.99 




b"- 




30.50llg:i 


29.88ilii:?i 


30.2011^:1 


31.4911^:°^ 


01 QQ+13.82 
Oi.OO_ii 48 


24.0 ±2.5 






32.00t}t?2 


Q1 10+14.60 


31.34li?:i 


00 00+14. 96 










n 47+0.20 




0.39l°:l 


0.301°:?^ 


0.301°:?^ 


1 oQ+0.56 
-'-•c>'J-0.55 








0.291°:! 


0.29l°:?f 


U-.^^-0.08 


n 90+0.20 

'-'•^'J-0.09 




B- - 






7.36111° 


7.47li:i 


7.ioli:i 


7 9-1 +3.21 
' •^J-_2.37 


6.6 ±0.6 








6.8411°^ 


6.90l|i 


6.54l|i 


6.60111 








U-Ui_o.01 


U-UZ„Q Ql 


09+0 03 
U-U^„0.01 


U.UUO„o.003 


004+°-0i^ 


< 1.2 






n 03+0-°^ 


09+0-02 


09+° °3 


010+°'°i'^ 

U.UlU_o.oo7 


010+° °20 

U.UiU_o.007 





with the experimental data within errors. Since these decay channels are dominated 
by the color-allowed tree amplitudes, both the weak annihilation and the higher order 
penguin contraction contributions are small. In addition, the main theoretical errors come 
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Table 4: CP-averaged branching ratios (in units of 10~^) of penguin-dominated (the 
upper six) and annihilation-dominated (the last two) B PV decays with AS = 0. The 
captions are the same as Table El 



Decay mode 



NF 



0{as) 0{as + a 



0{as 



a: 



EXP. 



B- 



'-'•-'-^-0.04 



0.181°:°? 



U-^<J-0.09 



< 5.3 



0.321°:}? 

0.141°:°^ 



9q+o-io 

U-'^'J-O.OS 

1 fi+°-°9 

U- J-O-0.06 





'J-'J'J-O.lO 

2fi+°-i2 



0.291°:}^ 

o.2ol°:J? 



0.401°:?° 



0.31 



^0.15 
-0.10 



B- 



B^ 



K^K*- 



*o 



0.29l°:J^ 

o.o6l°:J^ 
o.o5l°:Jl 
o.o6l°:J^ 



0.22 

o.io: 



-0.10 
-0.08 



h0.21 
-0.07 



0.08l°-i« 



-0.06 



0.091°-^° 



-0.06 



si+°-i^ 

0.071°:J^ 

o.o9l°:Si 



0.26: 

0.18: 



r^O.lO 
-0.09 



rO.27 
-0.10 



0.151°-^^ 



-0.09 



0.181°-^^ 



0.36; 

0.18; 



hO.16 
-0.11 

^0.26 
-0.10 



0.14l°-2' 



0.08 



-0.10 



0.17l°f 



-0.09 



04+°-^'' 



0.07; 



kO.16 
-0.05 



06+°-^^ 



1 5+°-25 

J-J-0.08 



14+0-24 
U--'-'±-0.08 



B- 



TX 



0.001 



0.008 



< 0.41 



TT 



i 0.001 
0.0003 
0.0003 



0.007 
0.004 
0.003 



< 1.0 



K*-K+ 



U.Ui8„o.004 



o.oi9l°:°g^ 



B^ 



K-K* 



U.Ui»„o.004 



0.019: 



hO.005 
-0.004 



from the uncertainties of the form factors and CKM matrix elements. 

• The decays B^ 7r°p° and B^ -k^uj. Since these decay channels are dominated by the 
color-suppressed tree amplitudes, their branching ratios are predicted to be very small. 
The higher order penguin contraction contributions are always much smaller than the 
weak annihilation contributions. Besides the form factors and CKM matrix elements, 
the spectator-scattering amplitudes also cause sizable uncertainties to their CP-averaged 
branching ratios. 
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Table 5: CP-averaged branching ratios (in units of 10 ^) of penguin-dominated B PV 
decays with AS* = 1. The captions are the same as Table El 



Decay mode NF 8^+" EXP. 









0{a,) 


0( 


as + 






9 \ 




B- - 


TT A 


„ in 70 

^•^'-0.64 


^•OU_o.88 


4.26l}i? 


f-n+i.22 


5, 


oq+2.01 


10.8 ±0.8 






4.89ii;i 


oc+i.^' 


5, 




, |_-| CI 


6, 


Q/i+2.18 
•<J^-1.70 




B- - 




-•-•0-^-0.54 


1 00 ~rU. 1 y 


2, 


70+1. 23 
• ' <J-0.81 


'^•'J'J-0.69 


3, 


•^^-0.89 


6.9 ± 2.3 








r) f) 1 i-U.o / 
^•^^--0.74 


3, 


nc-+1.25 
•UO-0.89 


^- ' ^-0.79 


3, 


70+1-36 
• ' U_l.oi 




B - 


n+K*- 


1.84l°|o 




3, 


•U'±-1.04 


^ J fc— Li OS 

2.47lJ:°^ 


3, 


70+1.84 
• '0-1.34 


- - — _Ll f\ 

ii.7tii 






, ^ _Ll /I Q 

3.40l|:n 


^•^^-0.84 


3.43t}:?I 


^ ^ ^ _Ll Q1 

2.991JJ1 


4.30^?:°^ 




B - 


TT^K 


0.49t[]:2^ 


U.OO„o 26 


1, 




^ ^ ^ -Lfl A 9 


1, 


4c;+0.86 
•^^-0.56 


1.7 ±0.8 










1, 


1 n 7Q 

OQ+U. ( 

•■^0-0.50 


i.07l°:i 


1, 


79+0.91 
• ' ^-0.65 




B- - 


K-(f) 


. _|_-| -I Q 

' -'--i.oo 


^' ' '^-1.20 


5, 


•UD-1.48 


. ^ < _l_1 f^Q 

4.04^1:1 


6, 


77+2.78 
• ' ' -1.74 


„ in 

9.03l°:i 






lo.iriti 


^ „ „ 1 -] no 

3.90^11^ 


6, 


00+2.07 
•'J^-1.77 


^ _L9 O"? 

5.59li:?? 


8, 


49+2.67 
•^^-2.22 




B - 


K (j) 


3.45tJi!i 


„ ^ 1 1 on 

2.53li:n 


4, 


7ri+1.86 
• ' U_i.37 


3.67tl:l? 


6, 


iq+2.40 


8.31lo 






9.46t|°^ 


„ |_-i 01 

3.63tl|} 


5, 


•^^-1.67 


_ „ ^ 19 1 n 

5.09l?i? 


7 7ri+2.55 
' • 'U_2.14 




B- - 


K p 




1 7/i+'i-uy 
' ^-1.16 


1, 


fi^+3-10 

•'J^-1.08 


„ 1 /I /1 9 

'3-iO-1.85 


3, 




< 48 






„ ^ „ _l_0 1 v 


1 Qfi+2.yy 

-'-•'JO-0.97 


1, 


90+2.69 
•^'-'-0.86 


2.731?:^ 


2, 


4q+3.83 

•^^-1.47 




B- - 






o.99lJ:5™ 


0, 


q«+i.69 
•yO-0.56 


-■-•^"-0.95 


1, 


CI +2.24 
•^-'--0.95 


4.23l°|? 






„ Li ^ 1 

0.58lJi^ 


o.78i^:i 


0, 


79+1.35 
• ' ^-0.36 


1 ^ n 

l.ZO_o.78 


1, 


.q+2.12 
•-'-^-0.70 




B - 






3.44lt:9? 


3 


01 +4.09 
•-J -"--LSI 


^•^ ' -2.67 


5, 


1 1 +5.18 
•-'--'--2.55 


^ ^ J_1 (\ 

9.9tl| 






^•^0-1.33 


3.04t?:^^ 


2, 


81 +3-77 


i ^ ^ J_f^ 1 Q 

4.86^11^ 


4, 


c,c,+5.00 
•33-2.32 




B - 


-T^ n 


i.42tj:?^ 


1 /^o+2 1 3 


1, 


Qn+2-12 

•yu_o.97 


/-,o+3 01 

3.03tt35 


2, 


•94t?:i 


5.1 ± 1.6 






1 09+1.79 
-L-'J^-o.re 




1, 


f;f5+1.97 


9 00+2.6I 


2, 


70+2.59 
• ' U-1.27 




B- - 








2, 


1 n+2.55 


3.071?:°^ 


2, 


qq+3.07 
.yy_i.44 


6.5 ±0.6 






0.40^°;?^ 


1 7P;+2.15 


1, 


«c:+2.27 
•"^-0.94 


2.6li?i° 


2, 


/1 7+3.25 
•^' -1.29 






K uo 


0.17t°:?^ 


-'-•U'j-o.es 


0, 


qq+1-67 


1 70+2.45 

-•-•'O-l.OO 


1, 


79+2.26 
• ' ^-0.96 


4.7 ±0.6 








n 7^+1.49 


0, 


.oy_o.47 


1 4Q+2.I6 
-'-•^'J-0.83 


1, 


QQ+2.09 
•'J'J-0.82 
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The decays K~ K*^ and B ^ K K*^. These decay channels are dominated by 

the h ^ d penguin amphtudes, and the dominant term is proportional to the coefficient 
a\{PV). Since 0^4 ~ and |A'„| ~ |A'^|, large interference effects between these two terms 
are expected and the branching ratios of these decay modes have a strong dependence 
on the weak phase angle 7 = arg{y*^). The higher order penguin contraction contribu- 
tions can provide about 60% enhancements to their branching ratios, and are larger than 
the weak annihilation contributions (which also play an important role in these decay 
channels). Since the higher order penguin contraction contributions are all involved the 
quantity A^, the main theoretical errors in the CP-averaged branching ratios, besides the 
CKM matrix elements and form factors, also originate from this quantity. 

The decays B^ K^K*~ and B^ — K'^K*^ . The dominant contribution to the decay 
amplitudes is proportional to the coefficient a\{yP), where delicate cancellations between 
the vector and scalar penguin contributions occur, their branching ratios are therefore 
predicted to be relatively small. This also renders the weak annihilation contributions 
potentially large. On the other hand, since these decay channels belong to the category of 
B VP decays, the higher order penguin contraction contributions are predicted to be 
small. The theoretical errors in the CP- averaged branching ratios of these decay channels 
are large, mainly due to the variations of the strange-quark mass and A^. 

The decays B^ 7i^(f) and P*^ n^cf). These two decay channels do not receive the weak 
annihilation contributions and are electro-weak penguin dominated processes. Due to the 
small coefficients a^{Ti(j)) and «3,e«;(^0); their branching ratios are predicted to be quite 
small. From the numerical results, we can see that large "nonfactorizable" contributions 
dominate these decays, while the theoretical predictions are still quite lower than the 
experimental upper bounds. The higher order penguin contraction contributions have 
negligible impact on these decay channels. 

The decays B^ K~^K*^ , K~ K*~^ . These two decay channels are pure annihilation 
processes. Studying on these decay modes may be helpful to learn more about the strength 
of annihilation contributions and to provide some useful information about final-state 
interactions. The higher order penguin contraction contributions have no impacts on 
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these decay channels. 

For penguin-dominated AS* = 1 decays, since the QCD penguin coefficients can receive 
large "nonfactorizable" contributions within the QCDF formalism, the predicted branching 
ratios for these decay modes are usually quite different from those obtained with the NF ap- 
proximation. In addition, the weak annihilation contributions to these decay channels are quite 
sizable. From the numerical results given in Table El we have the following general remarks: 

• The decays B ttK* and B K(j). With central values of our input parameters, our 
results are still lower than the experimental data. The dominant contribution to the decay 
amplitudes is proportional to the coefficient a4^{PV). The higher order penguin contrac- 
tion contributions can give enhancements to these branching ratios by about 40% ~ 90%, 
and reduce the discrepancies between the theoretical predictions and the experimental 
data. In addition, large interference effects between the tree and penguin amplitudes in 
some decay channels, such as 5^ 7r^K*~ and B^ n^K*^ , are expected. It is thus 
possible to gain some information on the weak angle 7 from these decay channels. The 
main theoretical errors are due to the uncertainties of the CKM matrix elements, form 
factors, and A^- 

• The decays B Kp and B Kuj. In their decay amplitudes, the dominant term is pro- 
portional to the coefficient aliVP). Because of the destructive interference between the 
vector and the scalar penguin contributions, the coefficient a\{yP) is reduced to a large 
extent, making the branching ratios of these decay modes much smaller than those of the 
corresponding B PP counterparts. It also makes the sub-leading terms, for example, 
the weak annihilation contributions, very important to account for the experimental data. 
Since these decay channels also belong to the category oi B ^ VP decays, the higher 
order penguin contraction contributions are quite small, and tend to decrease the NLO 
results. The main theoretical errors are due to the uncertainties of the strange-quark 
mass and form factors. 

From the above discussions, we can see that the higher order penguin contractions of 
spectator-scattering amplitudes play an important role in penguin-dominated B PV de- 
cays, while for tree-dominated B PV decays, their effects are generally quite small. In 
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particular, for decay modes dominated by the coefficient a^^lPV), these higher order penguin 
contraction contributions can increase the branching ratios by about 40% ~ 90%, while for 
those dominated by the coefficient a^{VP), their contributions are also predicted to be small 
and tend to decrease the branching ratios of these decay modes. At present, all these predicted 
CP-averaged branching ratios still suffer from large theoretical uncertainties. 

3.3 Direct CP-violating asymmetries of S ^ PV decays 

In this subsection, we will discuss the direct CP-violating asymmetries. In particular, we 
will investigate the impact of the higher order penguin contractions of spectator-scattering 
amplitudes on this quantity. 

We adopt the convention for direct CP asymmetries 

B(^^f)^B(B^ (42) 

30" ^ f) + B{B^ ^ f) 

Our numerical results for the direct CP-violating asymmetries are listed in Tables IHl and [7| Since 
the strong phases are suppressed by and/or Aqco/'^fc within the QCDF formalism, the direct 
CP- violating asymmetries for most B — * PV decays are predicted to be typically small. This is 
particularly true for decay modes dominated by the tree coefficient ai, for example, the decay 
P° 7r"p+. However, for h ^ d penguin dominated B — s> KK* decays, the penguin amplitudes 
generated by the internal u-quark loop and c-quark loop are proportional to the comparable 
CKM elements V*i,Vud and V*fjVcd, respectively, large direct CP-violating asymmetries for these 
decay channels are predicted. In addition, due to large interference effects between the tree 
and penguin amplitudes, the direct CP-violating asymmetry of B~ ti^K*~ decay is also 
predicted to be large. 

Due to cancellations among the strong phases associated with the individual Feynman dia- 
gram in Fig.inias discussed in Sec. 13. 11 the higher order penguin contraction contributions to the 
direct CP- violating asymmetries for most B PV decays are predicted to be small, however 
with a few exceptions. From Table El we can see that both the higher order penguin contraction 
and the weak annihilation contributions have significant impacts on the direct CP- violating 
asymmetry of P^ n^p^ decay. This is due to the delicate cancellations among the competing 
terms in its decay amplitude, making these sub-leading contributing terms potentially large. 
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Table 6: Direct CP-violating asymmetries (in units of 10~^) for two-body hadronic 
charmless B PV decays with AS* = 0. Decay modes with very small branching ratios 
are not considered. A(jp and denote the results without and with the annihilation 

contributions, respectively. The other captions are the same as Table El 
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From the numerical results given in Tables El and [3 we can also see that the higher order penguin 
contraction contributions to the direct CP-violating asymmetries of P^ ^r+p^, P~ tt'uj, 
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Table 7: Direct CP-violating asymmetries (in units of 10 ^) for two-body hadronic 
charmless B PV decays with AS* = 1. The captions are the same as Table El 
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B 7r^p+, and B decays are also quite large. In particular, these higher order pen- 

guin contraction contributions can increase the direct CP- violating asymmetries of the former 
two, while decrease those of the latter two by the same magnitude. 

Although the uncertainties from various input parameters are reduced to some extent, the 
renormalization scale dependence of the direct CP- violating asymmetries for some decay modes, 
such as B~ K~uj and B~ — > K~p^ decays, are still large. This is due to the fact that the 
imaginary parts of the coefficients define by Eq. (I10|) . which are crucial for the direct CP- 
violating asymmetries, generally have a larger scale dependence jH]. 

The direct CP- violating asymmetries for some hadronic charmless B PV decays have 
been measured recently, the data are still too uncertain to draw any meaningful conclusions 
from the comparison with the theoretical predictions, which also suffer from large uncertainties. 
With theoretical progresses and the rapid accumulation of experimental data, the situation will 
be improved and large direct CP- violating asymmetries in some decay channels, for example, 
the decays B KK*, will be found in the near future. 

3.4 Detailed analysis of B ^ ttK*, Kp decays 

The B ttK* and B Kp decays, like their PP counterparts B ttK decays, are also 
penguin-dominated decay modes, and hence sensitive to any new physics contributions. If the 
"tt/T" puzzles, with the improvement of experimental measurements, still remain unexplained 
within the SM, there would be signals of new physics beyond the SM j2ZI. Thus, the B ttK* 
and B Kp decays can be used to determine whether there are any new physics contributions 
and, which one, if exist as hinted by the "vrJi'" puzzles, is more favored. On the other hand, 
once the data on these decay modes becomes more precise, useful information on the weak 
phase angle 7 can also be obtained from these decay modes j2H] • So, detailed studies on these 
decay modes are worthy. 

In Figs, iniand^l we show the dependence of the CP-averaged branching ratios of these 
decay modes on the weak phase 7. In these two and the following figures, the central values of 
all input parameters except for the CKM angle 7 are defaulted and the renormalization scale 
is fixed at fi = nih. 

From these two figures, we can see that the experimental data on these decay modes are 
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Figure 9: The 7 dependence of the CP-averaged branching ratios of B ^ irK* decays. The upper and 
the lower four plots denote the results without and with the annihilation contributions, respectively. 
The solid and dashed lines correspond to the theoretical predictions with and without the higher order 
penguin contraction contributions, respectively. The horizontal solid lines denote the experimental 
data as given in Table 01 with the thicker ones being its central values and the thinner its error bars. 
The NF results denoted by the dash-dotted lines are also shown for comparison. 



generally larger than the theoretical predictions obtained based on the QCDF approach. Some 
decay modes, such as tt^K*~ and b'^ ir'^K*' decays, have a strong dependence on the 

weak angle 7. Moreover, both the higher order penguin contraction and the weak annihilation 
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Figure 10: The same as Fig. IHlbut for B Kp decays. 



contributions can give significant enliancements to tlie CP- averaged branching ratios of i? — >• 
irK* decays. However, for B tiK* decays, only the weak annihilation contributions can 
provide large enhancements to the CP- averaged branching ratios, and the higher order penguin 
contraction contributions play only a minor role. 

Since the theoretical uncertainties in the predicted CP-averaged branching ratios can be 
largely reduced by taking ratios among them, we shall discuss below certain ratios among the 
CP-averaged branching fractions of these decay modes, like the ones defined for B irK 
decays [22] ■ 
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For B ttK* decays, we define the following three ratios 



R{ttK ) = — — — 0- , (43) 

U„^^K^) . '-^'ir^'V- (45) 

With 7r{K*) meson replaced by p{K) meson, we can get another three similar ratios for B Kp 
decays. These ratios should be more appropriate to derive information on the weak phase angle 
7, as well as the relative strength of tree and penguin contributions than branching ratios. 

Our numerical results and the current experimental data for these ratios are presented in 
Table IHl The 7 dependence of these ratios are displayed in Figs. ITT] and IT^ 

From these two figures and the numerical results given in Table El we can see that our 
theoretical predictions for most of these ratios are in agreement with the data, considering the 
large uncertainties in the experimental data. 

From the explicit expressions of the decay amplitudes for these decay modes as given, for 
example, in Ref. ^2], we can see that differences between the two ratios Rc and Rn for both 
TiK* and Kp modes arise mainly from the color-allowed electro- weak penguin coefficient af 
and the color-suppressed tree coefficient 02, which are both predicted to be small within the 
QCDF formalism. So, the two ratios Rc and Rn are expected to be approximately equal 
within the SM. However, due to delicate cancellations among various competing terms, these 
ratios are strongly affected by the sub-leading contributing terms. After including the weak 
annihilation contributions, the two ratios Rc and Rn tend to be approximately equal. The 
current experimental data, however, indicate that RniirK*) is somewhat larger that Rci^irK*), 
but with large errors in the former. Unfortunately, due to the insufficient data on the branching 
ratios of the Kp modes, direct experimental comparison between Rc{pK) and Rn{pK) is not 
feasible for the time being. Once the experimental "-Rc — Rn' comparison with the case of ttK* 
and Kp decays are available, we can determine whether our theoretical predictions based on 
the QCDF approach are correct. 

It is also noted that the patterns of these quantities remain nearly unaffected even with 
these higher order penguin contributions included, because the higher order penguin contraction 
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Table 8: Ratios among the CP-averaged branching fractions of i? — > 7TK*,Kp decays. 
Numbers shown in columns 3 and 4 correspond to the results obtained without the an- 
nihilation contributions, while those in columns 5 and 6 the ones with the annihilation 
contributions. The other captions are the same as in Table El 
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contributions to the decays in the same ratio are similar in nature, and hence eliminated. 

With refined measurements available in the forthcoming years, it would be very interesting to 
check whether the theoretical predictions for these ratios are consistent with the data. Moreover, 
studies on these B PV modes will help us to understand the "vri^" puzzles j2Zj. 

4 Conclusions 

In this paper, we have reexamined the hadronic charmless B PV (with P = {tt,K), and 
V = {p, K* ,uj,(j))) decays in the framework of the QCDF. We have taken into account the 
penguin contractions of spectator-scattering amplitudes induced by the b Dg*g* transitions, 
which are of order a^. The main conclusions of this paper are summarized as follows. 
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Figure 11: Ratios of the CP-averaged branching fractions for B — > -kK* decays defined by Eqs. 

(|45j) as functions of the weak phase 7. The left and the right plots denote the results without and with 

the annihilation contributions, respectively. The meaning of the other lines is the same as in Fig. 



1. For penguin-dominated B PV decays, predictions obtained based on the QCDF ap- 
proach are generally quite different from the ones obtained with the NF approximation 
due to large "nonfactorizable" effects on the penguin coefficients. Contrary to their PP 
counterparts, the PV modes usually have smaller penguin coefficients af, rendering the 
sub-dominant terms potentially large. For example, the weak annihilation contributions, 
though power suppressed by AQCD/''^fe; are very significant in these penguin-dominated 
decays. The higher order penguin contraction contributions can interfere significantly 
with the next-to-leading order results, and hence are also important for these penguin- 
dominated decay modes. In particular, for decay modes dominated by the coefficient 
a\{PV), the higher order penguin contraction contributions can increase the CP-averaged 
branching ratios by about 40% ~ 90%, while for those dominated by the coefficient 
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Figure 12: The same as Fig. El but for B Kp decays 

a\{yP), their contributions are predicted to be small and tend to decrease the branching 
ratios of these decay modes. 

2. For tree-dominated decays and the decays having only the penguin coefficients a^, afg^ 
or having only the weak annihilation contributions, the higher order penguin contraction 
contributions to the CP-averaged branching ratios are predicted to be quite small. 

3. Since the direct CP-violating asymmetries are proportional to the sin of strong phase, 
which is usually suppressed by as and/or Aqcd/"2;, within the QCDF formalism, most of 
the hadronic charmless B PV decays are predicted to have typically small direct CP- 
violating asymmetries. However, for those decay modes where there are large interference 
effects between various contributing terms in the decay amplitudes, such as P ^ KK* 
decays, large direct CP-violating asymmetries are predicted. 

4. Due to large cancellations among the strong phases associated with the individual Feyn- 
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man diagram in Fig. IHl the higher order penguin contraction contributions to the direct 
CP- violating asymmetries for most B PV decays are predicted to be small, however 
with a few exceptions. For example, we find that both the higher order penguin con- 
traction and the weak annihilation contributions have significant impacts on the direct 
CP-violating asymmetry of b'^ 7r°p° decay. In addition, the higher order penguin con- 
traction contributions to the direct CP- violating asymmetries of P" Tr+p", B~ 7i~uj, 
P° — i> n~ 1 and P° K~ decays are also quite large. 

5. With more accurate experimental measurements available in the forthcoming years, it 
would be very interesting to check whether the theoretical predictions for the ratios P, 
Pc, and Rn for both the nK* and Kp decay modes are consistent with the experimental 
data. In particular, the experimental Pc — P„ comparison with the case of vrPT* and pK 
decays are very crucial for our understandings of the "vrP'" puzzles. 

Although the theoretical results presented here still have large uncertainties, the penguin 
contractions of spectator-scattering amplitudes induced by the h Dg*g* transitions, which are 
of order a^, have been shown to be very important for two-body hadronic charmless B — > PV 
decays, particularly for those penguin-dominated ones. It is very interesting to note that the 
1-loop (a^) correction to the hard spectator scattering in the tree-dominated B txtx decays 
has been performed recently [201, which forms another part of the NNLO contribution to the 
QCD factorization formula for hadronic P-meson decays. Using the PQCD method, the NLO 
corrections have also been carried out for P ttti^tiK and pp decays very recently In 
addition, much progresses in SCET have also been made in the past two years [321 • With the 
steady progress in experimental measurements at BaBar and Belle, further systematic studies on 
these higher order contributions to the rare hadronic P-meson decays are therefore interesting 
and deserving. 
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APPENDIX A: ANALYTIC EXPRESSIONS FOR THE 
Ai FUNCTIONS 

In the NDR scheme, after performing the loop momentum integration, subtracting the regulator 
e using the MS scheme, and performing the Feynman parameter integrals, we get the analytic 
expressions for the Aj functions appearing in Eqs. (jHBj) and (jHTjl 

2 ri 4 
Az5 = 2 + ^ [Go(ri)-Go(ri + r3)] -- [G_i(r0-G_i(ri + r3)] , (46) 

Ai, = _ 1 ^ 2n (1 + n) ^^^^^^ _ 2 (n + n.) ^^^^^ ^ 

+ -[G-M)-G-i{n + rs)] -^i^^^To(rO 

+ ''-''■-''■;"'•■ ^'•^' r„(n + r,), (47) 

2 ri 4 

A^23 = -2 i [Coin) - Coin + rg) ] + - [G_i(ri) - G„i(ri + rg) ] , (48) 

^3 ^^3 
A226 = -AZ23, (49) 

9 3 TTic 3r3 3r3 

+-[G_i(ri)-G„i(ri + r3)] , (50) 

^"3 



22 12 4ri 8, /i 2 (Tn - r3 - 3ri r3 + 2r2 - 2ri, 

Az3 = 77 + — + ^ gin Go{ri + rs) 

9 ^3 3r3 3 TTic 3r3 



r„(r.) - ^ <^ - - f <'•' + ro(r. + .3) . (51) 

^3 '^S 

32 16, /i 8(2 + ri)^, , 8 (2 + n + r3) ^ , , 
A.S = Y-yl-£;-^^Go(rO + ^^^Go(n+r3), (52) 



32 ^ 12 ^ 4ri ^ 16 ,^ /i ^ 2ri (7 + 2ri + 6r3) 
r 

2 (2r2-r3 (l-4r3)+ri (7 + 6r3)) 



Aii2 = -— + — + — + — In ^ + ^ -Coin) 

9 r3 3r3 3 3r3 



2 Go(ri + r3) 

[G_i(ri) - G_i(n + r3) ] + ^ '."^^^ Un) 

To To 
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r 



3 



-- [G_,(ri)-G_i(ri+r3)] , (54) 

^^3 



2 16 8n 2n{Arl + 3rs {8 + r-i)+n {20 + 7 rs)) 

Q ^ ^ 5~2~? ^ \ 1+^3 

3 rs 3r3 3r| (ri + ra) 



3r^ 



4 (4 - ri) ri 4 (4 - ri - rg) (ri + rg) 
2 ^o(n) + 2 ^ol'^i + '^s), (55) 

where the notations ri = /c^/m^, r2 = p^/m^, and = 2{k ■ p)/ml have been introduced. 
With mc replaced by m^, we can get the results for the b-quark loops. For light u, d, s quark 
propagating in the fermion loops, these Az functions can be evaluated straightforwardly. Here 
only the relevant Ai functions are given. Explicit expressions for the remaining ones can be 
obtained similarly. 

The functions Gi{t) and Tj(t) are defined, respectively, by 

Gi{t) = [^dxx'\n\l-x{l-x)t-i6] , (56) 
Jo 



X 



with the explicit form for To(t) given by jTH] 



4 arctan V— . 0<t<4 



To{t) = { ^ ^^P^ V ^ , (58) 

' 2i 7r+21n(v^-Vf^)-21n(v/f+v/t^) . t > 4 

while the explicit form for G_i^o(^) could be found in Ref. |17j . 



APPENDIX B: INPUT PARAMETERS 

In this appendix, we present the relevant input parameters used in our numerical calculations 
as follows. 

Wilson coefEcients. — The Wilson coefficients Gi{fi) have been reliably evaluated to next- 
to- leading logarithmic order fl^ I33j. Their numerical values in the NDR scheme at the scale 
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fi = mb {fih = V^hrrib ) are given by 

Ci = 1.080 (1.185), C2 = -0.180 (-0.367), C3 = 0.014 (0.028), 

= -0.035 (-0.062), C5 = 0.009 (0.011), Cq = -0.040 (-0.085), 
Cy/ae.m. = -0.009 (-0.029), Cg/ae.m. = 0.050 (0.107), Cg/a^.m. = -1.238 (-1.375), 
Cio/«c.m. = 0.243 (0.451), C^^ = -0.302 (-0.365), C^J = -0.145 (-0.169), (59) 



with the input parameters fixed as [21]: as{mz) = 0.1187, ac.m.(^w) = 1/129, m-yv = 
80.425 GeV, mz = 91.188 GeV, sin^ = 0.2312, mt = 172.7 GeV rub = 4.65 GeV, 
Ah = 0.5 GeV. 

The CKM matrix elements. — Here we use the Wolfenstein parametrization for the CKM 
matrix elements 



CKM 



1-^ 



A AX^{p - irj) 



V 



2 
i2 



-A 1 

AX^{l-p-ir]) -AX' 
and choose the four Wolfenstein parameters {A, A, p, and 77) as |2ZI 



AX' 
1 



0{X' 



(60) 



A — 0.825l;o;oj;g, 



A = 0.22622 ±0.00100, p = 0.207+^:[J^^, r/ = 0.340 ± 0.023, (61) 



with p and r] defined by p = p (1 — ^) , 77 = 77 (1 



Masses and lifetimes. — For the quark mass, there are two different classes appearing in the 
QCDF approach. One type is the pole quark mass which appears in the evaluation of the 
penguin loop corrections, and is denoted by mq. In this paper, we take 



= m^i = = 0, = 1.46 GeV, mh = 4.65 GeV. 



(62) 



The other one is the current quark mass which appears through the equations of motion and 
in the factor r^^. This kind of quark mass is scale dependent. Following Ref. ^2]; we hold 
{rriu + ^d){fJ')/^s{f^) fixed, and use rns{p) as an input parameter with the following values 



2m,(p)/(m„ + m<i)(p) = 24.2, m,(2 GeV) = (98 ± 20) MeV 

mb{mb) = 4.26 GeV |S3 , (63) 



where the difference between the u and d quark is not distinguished. 
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For the lifetimes and the masses of the B mesons, we choose 



TB^ = 1.643 ps [2E] 
tb^ = 1.527 ps 



mB, 



5279.0 MeV 
5279.4 MeV 



(64) 



as our default input values. The masses of the light mesons are also chosen from Ref. 

Light-cone distribution amplitudes (LCDAs) of mesons. — The LCDAs of mesons are also 
basic input parameters in this approach. In the heavy quark limit, the light-cone projectors 
for the B, the pseudoscalar, and the vector mesons in the momentum space can be expressed, 
respectively, as jHl 1121 

ifsmB 



Ml, 



(M|f)a/3 



4 



(1+^)75 {$f(0+ i- <^f(0} 



Q:/3 



4 
ifv 



(75 ^p{x) - /ip75i^^^ ^v^^) 



ki ■ k2 

fnv fv hh 



al3 



(65) 
(66) 
(67) 



fv ki-k2 J «/3 

where ki and k2 are the quark and anti-quark momenta of the meson constituents and defined, 
respectively, by 

P P 



^2 



(68) 



2xp-p' ' ^ ^ 2{l—x)p-p' 
It is understood that only after the factor ki ■ k2 in the denominator of Eqs. and ()67j] 
cancelled, can we take the collinear approximation, i.e., the momentum ki and k2 can be set 
to be xp and (1 — x) p, respectively, with p being the momentum of the meson. ^m{x) and 
^m{x) are the leading twist and twist-3 LCDAs of the meson M, respectively. Since the QCDF 
approach is based on the heavy quark assumption, to a very good approximation, we can use 
the asymptotic forms of the LCDAs 



$p(x) = $y(x) = 6x(l 



X) 



^pix) 



^^(x) = 3(2x- 1). 



(69) 



■^It should be noted ^1 that, in defining the light-cone projectors of light mesons, all three-particle 
contributions have been neglected. The leading-twist LCDAs are conventionally expanded in Gegenbauer poly- 



nomials 



$A/(a;,^) = 6a;(l - x) 



l + ^af(/.)Ci3/^)(2..-l) 



where the Gegenbauer moments are multiplicatively renormalized. The asymptotic form of the leading 

twist distribution amplitude is valid in the limit fi 00. With three-particle contributions being neglected, the 
twist-3 two-particle distribution amplitudes are then determined completely by the equations of motion. 
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With respect to the endpoint divergence associated with the momentum fraction integral over 
the meson LCDAs, following the treatment in Refs. [HI El]; we regulate the integral with an 
ad-hoc cut-off 



with Ah = 0.5 GeV. The possible complex phase associated with this integral has been ne- 
glected. 

As for the B meson wave functions, we need only consider the first inverse moment of the 



where the hadronic parameter Xb has been introduced to parameterize this integral. In this 
paper, we take Xb = (460 ± 110) MeV as our input value jH]. 

Decay constants and transition form factors. — The decay constants and the form factors 
are nonperturbative parameters and can be determined from experiments and/or theoretical 
estimations, such as lattice calculations, QCD sum rules, etc. For their definitions, we refer 
the readers to Refs. [111211110]. In this paper, we take the following numerical values for these 
input parameters 



= 130.7 MeV [M], /i^ = 159.8 MeV [M], /b = 216 MeV [12], 
fp = 205 MeV dS], U = 195 MeV [S], fx- = 217 MeV 
U = 231 MeV [IS], fp{l GeV) = 160 MeV [IS], /^(l GeV) = 145 MeV [IS], 
f^4l GeV) = 185 MeV [IS], /^(l GeV) = 200 MeV [IS], 
F_f^^(0) = (0.258 ± 0.031) [IS], F^^^'(O) = (0.331 ± 0.041) [IS], 
A^^^(O) = (0.303 ± 0.028) [IS], Al^^^\0) = (0.374 ± 0.034) [IS], 
A^^'^(O) = (0.281 ± 0.030) [IS], (72) 



where the form factors are evaluated at the maximal recoil region. The dependence of the form 
factors on the momentum-transfer can be found in Ref. [IS]. It should be noted that the 
transverse decay constant fy is scale dependent. 




(70) 



LCDA <I)f (0 defined by [S] 




(71) 
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APPENDIX C: THE GUAGE INDEPENDENCE 



In this appendix, we present a detail checking of gauge independence of our calculation. 
Firstly, we would check the gauge dependence of Fig. 5(b) with the gluon propagator 



1 

g2 



(73) 



where the factor —i6"'^ has been suppressed, and ^ is the gauge dependent parameter. 

Before the light-cone projectors for mesons are sandwiched, the scattering amplitude of this 
diagram is read as 



A oc 



Us{p->)(T^pqP{l +-f5)UbiPb) 



D"'^{e)D>''{q^), (74) 



where the spin indices, color indices, and SU (3) color matrices have been suppressed. It is easy 
to show that 



Vd{PlhaVd{P2) 



1 



(P2 -Pl)^ 
Vd{Pl)'l'^Vd{p2) 



Vd{Pl)l^Vd{p2) 



Vd{Pl){^2- ^l)Vd{p2) 
1 



/3 /3 ■ 
P2-P1 

{P2-P1Y 



(75) 



{.P2-P1Y' 

where k = p2 — Pi is the momentum of the gluon connected to the spectator d quark, pi and 
P2 are the momentum of d quark before and after scattering, respectively. In the last step, we 
have used the on-shell condition Vd{pi){i>i + rud) = {^2 + Tnd)vd{p2) = 0. It is also easy to 
show that 



Us{P5)(^f,pq''{l + IblUbiPb] 



D'-'iq') = \ 
q- 



Us{pb)(y''''qpO- + i5)ub{pb] 



Us{p->)(y^''q^,qp{l + l^)ub{pb) 



UsiP5)(^''^qpi^ + lb)Ub{Pb 



(76) 



From Eq. (|76|l . we can see that the gauge invariance of Qsg removes the gauge- dependent ^ 
term. 

From Eqs. (ff^ . (fTSj) . and (fTUj) . we can see that the scattering amplitude of Fig. 5(b) is 
independent of the gauge parameter ^. The above proof could be directly extended to that of 
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Fig. 5(c), Fig. 6(b) and (c), since the building blocks I^{k) and I^{k) defined by Eq|T7|are also 
gauge-invariant, i.e., k^I^{k) = k'^I^{k) = 0. 
For Fig. 5 (a), its amplitude reads 



A oc 



Vd{Pl)laVd{P2) 



Us{p5)(^fipq^{'^ + lb)Ub{Pb) 



X D""' D'"'' (p2) D^^'' (g2) V^,,.^. (g, p, fc) , 



(77) 



where V^va' (o', P, k) is the triple-gluon vertex. One can observe that the amplitude is gauge 
independent, because of 



ik"" 
ip" 



Vd{Pl)laVd{P2) 
Uu{Pz)lvVu{Pi) 



0, ik=p2-pi), 

0, {P=P3+P4), 



Us{p5)(^^lpq^{^ + l5)Ub{Ph] 



0. 



(78) 
(79) 
(80) 



Similarly, we can find that the amplitude of Fig. 6 (a) is also gauge independent. Using Eqs. ()78j] 
and (fTHj) . one can find that the amplitudes of Fig. 6(d) and (e) are gauge independent. 

In summary, we have shown that the amplitudes of the Feynman diagrams in Figs. 5 and 6 are 
gauge independent. The gauge independence of this subset Feynman diagrams is guaranteed by 
the on-shell external quarks and the gauge-invariance of O^g and 4-quarks operation insertions 
I^{k) and I^X^)- There are many 0{al) Feynman diagrams belonging to the group in Fig.4 
but not shown there. For example, the processes b sg* followed by g* qiQi ^ g* ^ uu 
or a gluon loop, which are gauge dependent separately. To keep gauge-independent, one must 
calculate the full set of Feynman diagrams in the category of Fig.4, which are not calculated in 
this paper. 
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